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New physics at the weak scale that can couple to quarks typically gives rise to 
unacceptably large flavor changing neutral currents. An attractive way to avoid this 
problem is to impose the principal of minimal flavor violation (MFV). Recently it 
was noted that in MFV only scalars with the same gauge quantum numbers as the 
standard model Higgs doublet or color octet scalars with the same weak quantum 
numbers as the Higgs doublet can couple to quarks. In this paper we compute the 
one-loop rate for production of a single color octet scalar through gluon fusion at the 
LHC, which can become greater than the tree level pair production rate for octet 
scalar masses around a TeV. We also calculate the precision electroweak constraint 
from Z — > bb] this constraint on color octet mass and Yukawa coupling affects the 
allowed range for single octet scalar production through gluon fusion. 
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I. INTRODUCTION 

The standard model for strong, weak, and electrodynamics has been very success- 
ful. Nonetheless, most physicists expect that the next generation of very high energy 
accelerator experiments that are able to produce resonances with masses around the 
TeV scale will find physics beyond what is in the minimal standard model. There 
are two reasons for this. Firstly, the mechanism for weak symmetry breaking in 
the minimal standard model, a single scalar doublet, is the simplest but has no di- 
rect confirmation from experiment. Secondly, an awkward fine-tuning of parameters 
must be made, order by order in perturbation theory, to keep the physical mass of 
the Higgs scalar very light compared to the ultraviolet cutoff. This is called the 
hierarchy problem and is the motivation for most of the proposed extensions of the 
minimal standard model with new physics at the weak scale. 

Models with new degrees of freedom at the weak scale typically give rise to 
unacceptably large flavor changing neutral currents (FCNC) if the new physics 
can couple at tree level to the quarks and the most general couplings are al- 
lowed. The minimal standard model does not have such a problem because of the 
GIM mechanism. More generally, the large FCNC problem does not arise if the 
SU{3)q^ X SU{3)uji X SU{3)Dji quark flavor symmetry is only broken by a single 
pair of Higgs Yukawa coupling matrices g'^ and g^. This way of suppressing FCNC 
is called minimal flavor violation (MFV) In MFV, only scalars with the same 
gauge quantum numbers as the Higgs doublet or color octet scalars with the same 
weak quantum numbers of the Higgs doublet can Yukawa couple to the quarks [2]. 
Therefore it is of interest to study the properties of such scalars. Models with sev- 
eral color singlet, weak doublet scalars have been extensively studied. In [3|], certain 
decays of color octet scalars were studied in the context of Pati-Salam unification. 
Here, we continue the phenomenological analysis of color octet scalars begun in 
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The tree level pair production cross section for charged or neutral color octet 
scalars was computed in [2]. The precision electroweak variables S and T were also 
computed^ and a number of FCNC processes were considered. Here, we calculate 
the one-loop production rate for a single neutral octet scalar through gluon fusion, 
gg S^. We also derive the constraint on the strength of the coupling of color octet 
scalars to up-type quarks that arises from experimental data on Rb. This constraint 
restricts the magnitude of the the gg —>■ cross section. 



II. THE MODEL 



The standard model quark Yukawa couplings of the quarks to the Higgs doublet 
H are 

L = -g^jURiQjH - g^jdmQjH^ + h.c., (1) 

where i and j are flavor indices, and gauge indices have been omitted. Repeated 
flavor indices are summed over. The Yukawa couplings generate the mass matrices 

MS=9Uh'), MS = gS{H^y, (2) 

for the charge 2/3 and —1/3 quarks when the Higgs field gets a vacuum expectation 
value {H^) = 0, {H"^) = v/^/2. In the minimal standard model the only couplings 
that violate the SU (3)q^ x SU (3)(7^ x SU (3) Dr quark flavor symmetry are the Yukawa 
matrices g^ and g^. We can view the theory as being invariant under the flavor group 
if the Yukwawa matrices are endowed with the transformation property, 

g"" - Vug'^V^ g"" - Vng^'V^, (3) 



^ The parameter U is small. 
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where Vu is an element of SU{3)uji, Vd is an element of SU{3)Dji and Vq is an 
element of SU{3)qj^. 

In this paper we add to the minimal standard model a single weak doublet of 
color octet scalars S^. According to MFV, its Yukawa couplings to the quarks are, 

L = -Vu97j^R^T''QjS^ - VDg?,dR{r^Q^S^^ + h.c. , (4) 
where the matrices and also transform as 

f - Vvfvl - Vng^'vl (5) 

and are composed from g^ and g^ . So, 

f = 9'' + e^^9''{9'')'9'' + ... (6) 

and 

Note that there is a term in the ellipses of Eq. ([6]) proportional to g^ig^Yg'^] however, 
it does not give rise to flavor changing neutral current effects, so we neglect it. We 
will assume that the e'^'^ are small and that terms with more powers of the Yukawa 
couplings are more suppressed and can be neglected. Diagonalizing the quark mass 
matrices we find, in the quark mass eigenstate basis, that the couplings of the octet 
scalars take the form. 



- ^rjndm'^{s,,+2e?Vi{m'j!/vyVk, + ... ^T^dLjS^'^ 

- x/2r/z)rfffi— (l + 2ef(my/t;)2 + ... ) \/.J.TV,5^- + h.c 
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where V is the CKM matrix and mf , mf are the charge 2/3 and charge —1/3 quark 
masses. Since the top quark is by far the heaviest quark we can approximate the 

above by, 

V V ■' V 

- ^nndR^ (5,, + 2e^Vl{mt/vfV^, + ... ) T^'d^jS^'^ + h.c. (9) 

The term proportional to ef gives rise to a tree level contribution to flavor changing 
neutral current processes like BB mixing from 5" exhange. However, the leading 
contribution to B -^57 does not involve ef and experimental data on this process 
provides an important constraint on rjD. The parameter rjjj is constrained from data 
on the precision electroweak variable Rb] this is discussed in the next section. 
The most general renormalizable scalar potential is [2I, 

y = - ( H^'Hi - y ) + 2m|Tr S^'S^ + \iH^'H.{Ii S^^S^ + X2H^'HjTi S^^ Si 

+ h^H^'H^^Ti SiSj + X^H^'Tt S^^SjSi + X^H^'Ti S^^S.Sj + h.c. 
+A6Tr S^'SiS^^Sj + ArTr S^'SjS^^S^ + AgTr S^'SiTi S^^Sj + AgTr S^'SjTi S^^ Si 
+AioTr SiSjS^'S^^ + AnTr SiSjS^^S^\ 

(10) 

We have explicitly displayed the SU{2) indices on the Higgs doublet and on the color 
octet scalars. Traces are over color indices and the notation S = S^T^ is used, where 
the SU{?>) generators have their standard normalization: Tr T'^T'^ = jl. The 
coupling A3 has been made real by a phase rotation of the S fields. With this phase 
convention the phases of riu^D and A4,5 represent additional sources of CP violation 
beyond those in the minimal standard model. 



6 



The Higgs vacuum expectation value causes a tree level mass splitting between 
the octet scalars. It is convenient to decompose the neutral complex octet scalars 
into two real scalars, 



5- = ■ (11) 



Then the tree level mass spectrum is 



V 



2 
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rrigo = m| + (Ai + A2 + 2A3; ^ , 



,2 



mi.0 = m's + (Ai + A2 - 2A3) —. (12) 

In this paper we focus on color octet scalars with masses greater than 500 GeV. 
The mass splittings are expected to be small compared with this and so we neglect 
the mass splittings between the various color octet scalar states for the remainder of 
this paper. Color octet scalars with masses between SOOGeV and ITeV can have a 



dramatic impact on the rate for Higgs production at the LHC 



a. 



III. PRECISION ELECTROWEAK CONSTRAINTS 

In [2j the values of the oblique parameters S and T that arise in this model 
were computed. These corrections to the standard model are expressed in terms of 
parameters in the scalar potential (including the color octet scalar masses). The 
contribution to the effective Hamiltonian for 6 ^ 57 proportional to rjurjo was also 
considered. Here we complete the analysis of precision electroweak physics in this 
theory by computing Rh, the ratio of the Z width to final hadronic states containing 
a b and b quark to the total hadronic width. We write the coupling of the Z boson 
to quarks q as, 

Z>^q^,{h,,PL + fR,,PR)q, (13) 



cos 6* 



w 



7 



b b 




FIG. 1: Feynman diagrams contributing to Zbb vertex correction. 

where Pl,r are the projectors Pl = {1 — 75)/2 and Pr = (1 + 75)/2. We use a 
supercript "0" to denote the tree level standard model value of the coupling, 

fl^ = tl~ sin'OwQ, fl, = -sm'9wQr (14) 

We take the value of siia^Ow from the measured vector and axial vector lepton cou- 
plings at the Z pole; then corrections to Z vacuum polarization effects from the octet 
scalars are absorbed into it. 

The down Yukawa coupling parameter, 7]^, is constrained by the B Xg'y partial 
width to be significantly less than mt/nify when r/^ ~ 1 and ms ^ 5TeV. Therefore, 
we neglect the down Yukawa coupling to the b quark, as well as neglecting quark 
masses other than the top. Then the octet scalars give a one-loop correction to 
the left-handed bottom coupling from the vertex diagrams and bottom quark wave 
function renormalization diagrams pictured in Fig. [H Writing, 

fL,b = + Sh,h, (15) 
it is conventient to decompose the correction to the left-handed bottom quark cou- 
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pling as, 



We find that, 



1 



3 V 167r2 



(16) 



pi pl—x 

A = 2 dx dy 
Jo Jo 

+ 



-In 



m^x + — x) 



m|(x + y) + m^{l — x — y) — M^xy 



mr 



ml[x + y) + TTigil ~ X — y) — M^xy 



(17) 



and 



pi p1—x 

B = 1 + 2 dx dy In 
Jo Jo 

ml/2 + xyMl 



m1{x + y) + — X — y) — M'^xy\ 
m|(x + y) + rn1{l — x — y) — M^xy J 



m1{x + y) + w^|(l — X — y) — M^xy 
Here ms is the mass of the charged octet scalars. Note that in the hmit ms » 

mt» Mz, 

2 / 2 \ 

(19) 



A = 2B = 24ln(4\- 



ml. 



m 



We compute 6Rb using the formula, 



L,b 



-0.78 5fL,b 



(20) 



where 5fL,b comes from the vertex and h quark field renormalization diagrams in 
Fig. [T] and is given above. 

According to the Particle Data Group , the observed value for at the Z pole 

is 

Rb = 0.21629 ± 0.00066 (experiment) (21) 
and the Standard Model predicted value is 



Rb = 0.21578 ± 0.00010 (SM prediction). 



f22l 
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FIG. 2: One (solid line) and two (dashed line) standard deviation exclusion contours due 
to Rfy. (Parameter space above these lines is excluded.) The curves were calculated using 
mt = 170.9 GeV, sin^ 9^^, lept = 0.23153, mz = 91.1876 GeV, and v = 246 GeV. 



Noting that the right-hand side of Eq. ( ]T6il is positive, the following one sigma bound, 

Sf^ < 0.00020, (23) 

on the contribution of the octet scalars to 6f^ follows from comparing the experi- 
mental value and the standard model prediction for Rb. 

In Fig. [2] we plot the one and two sigma curves that bound the excluded Irjul-ms 
parameter space. For example, if rris = ITeV, then \7]u\ must be less than 0.8 to 
agree with experiment to within la, or less than 1.8 to agree to within 2a. 



IV. THE CROSS SECTION FOR SINGLE COLOR OCTET SCALAR 
PRODUCTION VIA GLUON FUSION 

There are two neutral color octet scalars in the model we are considering. They 
are the scalar and pseudoscalar states destroyed by the real and imaginary parts of 
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the field S^^. The one loop gluon fusion rate p for producing these states singly is 
related to their two gluon decay rates by the standard formula,^ 

a{pp^S'},^jX)=T{Slj^gg)d^) -g{x)g{my sx), (24) 

\sms J Jjnl/s X 

where ^ = 1/16 is a spin-color factor that takes into account the interchange of 
summed-over and averaged-over states in the production cross section and decay 
rate. 

We begin by discussing the color octet scalar state's production. To simplify the 
analysis we assume that CP nonconservation is small and take t]u, A4^5 to be real. 
Neglecting the mass differences between the various charged and neutral color octet 
scalar states we find that 



nS'n - 99) = 



V2 2-7r 



10^3 



3C,M^, + A5)^(y - l)Re/(m?/m|) + ^dX, + A5)^^(y - l)^ 



■ (25) 



In Eq. (1251) . I{z) is the familiar factor from standard model Higgs decay. Assuming 
2; < 1/4 it is given by, 

2 

/(.) =2z + z{4z - 1)^ f{z) = (in (}^^^) - • (26) 
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The factors Cj in Eq. (125|) are the color factors, 

C, = Y^id^^^^f = C2 = Y1 d^'^'^d''^'^ f^^^ /^^^ = -20, (27) 



^ See, for example, 

3 The unique symmetric invariant, d,-^^^ , ofSU{n) for n > 3 is given by d^^^ = 2Tr ({T^, T^jT*^) 
where T"^ are the fundamental representation matrices. Recall that the structure constants, 
fABC^ are given by /^^^ = -2^Tr ([T^^, r^]T^). 
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FIG. 3: Diagrams contributing to 5"^ production via gluon fusion. For real A4^5, only the 
top loop contributes to production. 

and 

C3 = fAEFfBGEy ^ ^28) 

The last two terms in ( l25l) . which come from the octet scalar loops, are much smaller^ 
than the top loop contribution. (See Fig. O) This is partly due to the factor of 
(tt^/Q — 1). For the pseudoscalar we find, 

r(5y ^ gg) = ^^-——Cir]u \f{mjms)\ , (29) 

where f{z) is as in (126|) . The pseudoscalar rate is due solely to a top loop and is 
related to that of a heavy color singlet by a simple multiplicative factor. 

Allowing CP violation mixes the scalar and pseudoscalar states. For example, 
the induced ggS'}^ coupling acquires an axial contribution (and the ggSj coupling a 
non-axial contribution) proportional to lm{riu). Also, the 5*^ decay rate picks up a 
neutral scalar loop contribution proportional to Im(A4 + A5). 

In Fig. m we plot the cross section for single and 5*^ production, and real 
neutral scalar pair production as a function of m5, using the above results and the 
pair production result from [2]. For this plot the X^^^ terms are neglected and rju is 
set equal to unity. We expect that, just like in standard model Higgs production. 



^ This is true for A4.5 ^ 1. 
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FIG. 4: Production cross sections (in femto-barns) at LHC center of mass energy y/s = 14 
TeV for two real neutral scalars (solid line: pp S'^^j^S'^^j^X), one real neutral scalar 
(long dash: pp — > S^X), and one real neutral pseudoscalar (short red dash: pp — > S^X). 
The single color octet production cross sections are plotted with r]ij = 1 and include only 



igible for A4^5 ~ 1. We 
, and we used the two- 



the top loop contribution, as the scalar loop contributions are neg 
used CTEQ5 next-to-leading order parton distribution functions [8 
loop f3 function to run as{mz) = 0.1216 up to 0^(27715') for scalar pair production and 
as{ms) for single scalar production. The curves were calculated using mt = 170.9 GeV. 

the higher order QCD corrections are significant j^, [l^, |ll|. At a TeV, the single 
production rates begin to dominate over pair production. 

In Fig. we plot the one-loop single production cross sections for S'^ and Sj at 
the Tevatron. The tree level cross section for pair production is more than an order 
of magnitude smaller than the single production values in this mass and energy 
regime. For example, at ms = 500 GeV, a{pp — > S'^f^j^S^^j^X) = 10~^ fb, while 
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FIG. 5: Production cross sections (in femto-barns) at Tevaton center of mass energy 
^/s = 1.96 TeV for one real neutral scalar (solid line: pp Sj^X), and one real neutral 
pseudoscalar (dashed line: pp — > SjX). For this plot, r]u = 1 and the A4^5 terms in (|25|) are 
neglected. We used CTEQ5 next-to-leading order parton distribution functions [8], and 
we used the two- loop /? function to run as{mz) = 0.1216 up to Usims)- The curves were 
calculated using rrit = 170.9 GeV. 

a{pp Six) = 0.5 fb and a{pp -> S^X) = 0.9 fb. 



In a recent paper 121], the CDF collaboration presented preliminary limits from 
Run II on the production cross section for a "Z'-like" heavy neutral boson times 
its branching ratio to ti pairs. The analysis assumes that the boson appears as a 
Lorentzian enhancement in a limited region of the M^f spectrum. Their preliminary 
limits do not exclude neutral octet scalars that decay mostly to tt. 
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V. CONCLUDING REMARKS 

Minimal flavor violation is a convenient way to suppress flavor changing neutral 
currents when new degrees of freedom at the weak scale couple to quarks. If there 
are new scalar resonances with masses at the TeV scale that couple to quarks, then 
minimal flavor violation implies that they are either color singlets or color octets 
with the same weak quantum numbers as the Higgs doublet. Models with two (or 
more) Higgs doublets have been studied extensively. The phenomenology of modes 
with an additional color octet scalar were studied in j^. Here we extended this work, 
calculating the constraint on the strength of the Higgs coupling to up-type quarks 
that arises from precision electroweak data on We also computed the rate for 
single octet scalar production through gluon fusion. For color octet scalars with 
masses greater than 1 TeV this one-loop process can dominate over tree level pair 
production at ^/s = 14 TeV because the gluon parton distribution function increases 
rapidly as the momentum fraction decreases. 
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